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Abstract. Ever since we have been in the possession of quantum theories without observers, such as 
\ Bohmian mechanics or the Ghirardi-Rimini-Weber (GRW) theory of spontaneous wave function 

collapse, a major challenge in the foundations of quantum mechanics is to devise a relativistic 
quantum theory without observers. One of the difficulties is to reconcile nonlocality with relativity. 
I report about recent progress in this direction based on the GRW model: A relativistic version of the 

■ model has been devised for the case of noninteracting (but entangled) particles. A key ingredient 
\ was to focus not on the evolution of the wave function but rather on the evolution of the matter in 

V^J ■ space-time as determined by the wave function. 
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: INTRODUCTION 

o 

^ I I describe some models I designed recently that generalize the Ghirardi-Rimini-Weber 

(GRW) model [1, 2, 4] of spontaneous wave function collapse. One of these [5] is a 
relativistic version of the GRW model, the other [6], which I will describe first, an 
extension to quantum field theory (QFT). It would be nice to combine the two models 

■ into a single one suitable for relativistic quantum field theory; I have not started this task 
^ ! yet, but I hope I will be able to do this soon. In this article, I will begin by elucidating 
qh| the flash ontology for the GRW model [2], which is crucial for understanding the two 
^ new models. 
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: THE FLASH ONTOLOGY 

What are tables and chairs made of? To this question, different theories, even different 
versions of the GRW model, may propose different answers. Classical physics would 
say particles, described mathematically by their world lines. The objects that the answer 
describes, whatever it may be, were termed the "primitive ontology" of the theory by 
Durr, Goldstein, Zanghi, and AUori [7, 8, 9] and the "local beables" by BeU [10]. To 
give you concrete examples to think of, let me introduce two possible answers that have 
been proposed for the GRW model: 
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A) Not the wave function, but a continuous distribution of matter with density 

A' 

E 
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= j dVi---d3rr--dVAf ^/(ri, . . . , r?,r, . . . , ta^,? 



(1) 



where the hat denotes omission. This is (a simplified version of) a proposal due to 
Benatti, Ghirardi, and Grassi [11]. 

B) Not the wave function, but discrete "flashes", i.e., elementary events represented by 
space-time points. This is a proposal due to Bell [2, 12, 13]. In such a world, matter 

consists of lots of dots in space-time (see Fig. 1), and "a piece of matter then is a 
galaxy of such events" (Bell [2]). Such a space-time point Z = (R, T) is defined in 
the GRW model as the center of a wave function collapse. 
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FIGURE 1. A typical pattern of flashes in space-time, and thus a possible world according to the flash 
version of the GRW theory 

I will use the flash ontology in the following. It is misleading to say that the flashes 

are the collapse centers, though the collapse centers determine where in space-time the 
flashes are. It is better to think of the flashes as what the theory is about; when simulating 
the GRW model with the flash ontology on a computer, the set of flashes is the output 
of the algorithm. That is also why I prefer the word "flash" to "hitting," which would 
suggest instead that the wave function be the output of the theory. The right way to think 
of the primitive ontology is as describing the motion of matter in space-time, and thus 
as analogous to the particle trajectories in classical (or Bohmian) mechanics. A more 
thorough discussion of the notion of primitive ontology can be found in [9]. 

The flash ontology is an unusual choice of ontology; a more normal choice would be 
particle world lines, or fields in space-time. The motivation for this choice lies in the fact 
that the GRW model with flashes can be made Lorentz invariant (by suitable corrections 
in the equations), whereas the GRW model with the matter density m(r, t) cannot in any 
known way. The GRW model without primitive ontology, where only the wave function 
is regarded as existing, suffers, I think, from serious problems [9]. The basic problem is 
that a theory that never talks about matter would not form an adequate description of our 
world. Strictly speaking, in a world governed by such a theory there exists no matter, or 
else the link between matter and the mathematical variables of the theory is left vague. 
A related problem is that it remains to be clarified in what way the wave functions i/z^ 
and x^fu associated with different spacelike 3 -surfaces Z and YJ have to be compatible to 
be regarded as describing a consistent reality. In case L and H! have a large portion in 
common, one would presumably require that i/A^ and i/A^/ describe the same reality on 
ZnZ'. But the heart of this problem is that one has not defined clearly what that reality 



is and how the wave function influences that reality, and to provide such a definition is 
exactly the purpose of the concept of primitive ontology. 

When we regard GRW as a theory about flashes then the flashes form a random set, a 
point process (as mathematicians would say), in space-time. The joint distribution of the 
flashes is then determined by the (initial) wave function. Indeed, for this distribution one 
can write an explicit formula, which I will do in the next section, albeit in much greater 
generality including quantum field theories. 



GRW AND QFT 

It will be useful to write the abstract mathematical structure of the GRW model in the 
following way [6]. The mathematical ingredients are: 

• The initial wave function (at time to), a vector X]/ in a Hilbert space Jif with 
= 1. 



• For writing the flash rate (the probability per time with which a flash occurs) at a 
location r e in the form 

(V^lAWk), (2) 

we need a positive self-adjoint operator A(r) for every r e R^. 
For the case of nonrelativistic quantum mechanics of distinguishable particles, it 
is useful to introduce A'^ types of flashes (corresponding to the A'^ variables r,- G 
of the wave function), and thus for every type j e {1, . . . , A/^} the flash rate operators 
A,(r). To obtain the original GRW model, choose the flash rate operators to be 
multiplication operators by (3 -dimensional) Gaussians, 

A,(r) Vin ,...,r,)='- j^^, yfin ,...,r^), (3) 

where the model constants a and T are the width of the localization (say, a = 10^^ 
meters), and the average time between two collapses of the same type / (say, T = 10^ 
years). 

• The Hamiltonian H defines the way X]/ would evolve if there were no collapses. 

Together, J^^H, and A define a flash theory by the postulate that the joint distribution 
of the first n flashes is given by 

¥{Xi e d\, ...,x„e d\) = \\KnXi/fd\ ■ ■ -dV , (4) 

where d'^Xk = d^ri^^dtk is the 4- volume, Xk= {vkih), and 

Kn = A(r„)i/2 W,„_,„_, • • ■A{n)"^W^,.t, (5) 

with 

^ ^ fexp(-i/fr/^- i/d3rA(r)?) ift>0 
' 10 ifr<0. 



It is an easy exercise to check formally that with these definitions, (4) consistently 
defines a probability distribution on the suitable history space (a history in the flash 
ontology is a pattern of flashes, and history space thus is the space of all discrete subsets 
of space-time). The original GRW model is a special case of this scheme, for which 
/ d^r A(r) is, by (3), a multiple of the identity operator in so that Wt is the unitary 
Schrodinger evolution up to a scalar factor exp(— r/2T) that governs the (exponential) 
distribution of the random waiting time until the next flash. For the case of several types 
of flashes there exist obvious analogs to equations (4)-(6), see [6]. I also mention that 
equations (4)-(6) fit into a general mathematical scheme that Blanchard and Jadczyk 
[14] have developed for a different purpose, namely for models of coupled classical and 
quantum systems; see Jadczyk's contribution to this volume for details. The subsequent 
work [15, 16] on relativistic extensions of their approach, in contrast, seems not related 
to any relativistic theory of flashes. 

For A'^ identical particles, one should form the operators 

!=1 

restricted to the subspace of anti-symmetric (or, symmetric) wave functions for fermions 

(respectively, bosons). We thus obtain, in a remarkably simple way, a variant of the GRW 
model appropriate for identical particles, despite the suggestion sometimes conveyed 
in the literature (see, e.g., pages 312 and 382 in [4]) that for dealing with identical 
particles one has to employ a "CSL" model (continuous spontaneous localization [17]), 
in which the state vector follows a diffusion process in Hilbert space and thus "collapses 
continuously" rather than at discrete times. In fact, the modification of the GRW collapse 
process corresponding to the choice (7) of flash rate operators was already proposed in 
1995 by Dove and Squires [18], but has not received the attention it deserves. 
For a variable number of particles, one should form 

oo 

A(r) = 0AW(r) (8) 

on Fock space. In a sense, quantum field theory is nothing but the quantum mechanics 
of a variable number of identical particles. This allows us to define a flash process for 
more or less every QFT that is regularized (i.e., possesses a well-defined Hamiltonian). 
Equivalently, we can write A(r) as smeared-out "particle number density" operators 
N{r). 

(In nonrelativistic model QFTs, N{r) can be expressed in a simple way in terms of the 

field operators (j){r) of a quantum field, N{r) = 0*(r) 0(r). Since Pearle and Squires [19] 
have argued for choosing the collapse rate proportional to the mass, a better alternative 
to eq. (9) may be to replace the particle number density operators by the mass density 
operators.) 

This model is as hard to distinguish experimentally from standard QFT as the GRW 
model from standard quantum mechanics. Thus, contrary to what is sometimes conveyed 



in the literature, it is not a matter of necessity to employ a CSL model for dealing with 
QFT. I would rather say that the generalization to QFT of the GRW model with flashes 
is simple and natural. 

This model does not allow superluminal signalling by means of nonlocal correlations, 
and does not lead to any divergences. Ref. [20] includes a discussion of superselection 
rules in this model. 

Where is the collapsing wave function in this model? With every time t there is 
associated a wave function 

where n is the number of flashes between Iq and t, all of which enter into Kn. This 
function i/A; is what is best regarded as the wave function at time t; it is called the 
conditional wave function, since it depends on the realization of the flashes between tQ 
and L It is also called the collapsed wave function since each flash contributes a factor to 
Kn that collapses ^f (in one variable) to a neighborhood of the flash. As time proceeds, 
i/^f collapses whenever a flash occurs, and evolves deterministically in between. As a 
consequence of (4), the distribution of the flashes after t, conditional on the flashes 
between to and t, is given by the formula (4) for P with ^ replaced by % and ?o replaced 
hyt. 

Another fact I would like to mention for later reference is that the probability that no 
flash occurs before time t is given by 

Wt-tM?. (11) 



RELATIVITY 

A relativistic quantum theory without observers must overcome the difficulty of recon- 
ciling relativity with the quantum nonlocality discovered by Bell [3], which asserts that 
sometimes spacelike-separated events influence each other. 

A radical explanation of nonlocality postulates the existence of a preferred slicing, 
or foliation (as mathematicians would say), of space-time into spacelike 3 -surfaces. 
This foliation might naturally arise from the space-time geometry; for example, it might 
consist of the 3-surfaces of constant timelike distance from the big bang. Assuming a 
preferred foliation, a simple and natural extension of Bohmian mechanics to relativistic 
space-time has been worked out by Diirr et al. [21], and something analogous can 
presumably be done for the GRW model [22], be it with the matter density or with 
the flash ontology, see Fig. 2. 

While the existence of a preferred foliation seems a possibility worth considering, 
the conventional understanding of relativity forbids it. It is hard to make precise what 
is allowed and what is not, but it would seem that special 3-surfaces playing the role 
of simultaneity-at-a-distance are not allowed. If a preferred foliation exists, the conven- 
tional understanding of relativity is wrong. This motivates the search for alternatives, 
i.e., for relativistic theories that do not require special foliations or any other additional 
structure in space-time providing a synchronization of spacelike- separated events. Such 
an alternative is not in sight with Bohmian mechanics, but such a version of GRW with 




FIGURE 2. Given a foliation of space-time into spacelike 3-surfaces, one can define a version of 
the GRW model with matter density ontology, in which a collapse of the wave function leads to a 
discontinuous change of the matter density m{r,t) (shades of grey) along a 3-surface Z from the foliation. 



flashes is possible! 1 will describe one in the next section. This confirms a hope expressed 
by Bell in 1987 [2]: 

1 am particularly struck by the fact that the [GRW] model is as Lorentz 
invariant as it could be in the non-relativistic version. It takes away the ground 
of my fear that any exact formulation of quantum mechanics must conflict 
with fundamental Lorentz invariance. 

The observation that induced Bell to write this was the multi-time translation invari- 
ance of the original GRW model: Consider two quantum systems at a great distance 
from each other (far greater than their sizes); suppose the two systems are entangled but 
do not interact, H = Hi + H2, and that each is adequately described by nonrelativistic 
quantum mechanics. Then the effect of a small Lorentz boost is, to first order, a relative 
time translation; that is, the two systems get shifted in time by different amounts. In other 
words, exactly the synchronization (of events in one system with those in the other) that 
a preferred foliation would provide gets changed. The claim is that the joint distribution 
of the flashes is covariant under relative time translations, provided the two systems have 
different types of flashes. To see this, one may start from the following formula (here it 
suffices to consider = 2): It is the analog of formula (4) for noninteracting systems 
with different types of flashes, and provides the joint distribution of the first ni,. ..,nN 
flashes with label l,...,N. We denote by Xji^ the k-th flash of type i. 



where the operator A^, ,„ is given by (5) and (6) with A replaced by A, and H replaced 
by Hi. Now consider a pattern of flashes and apply a time translation to system 1 alone, 
shifting all flashes of type 1 to the past by some amount A > but leaving the flashes of 
type 2 unchanged. Apply this shift to the probability distribution and conditionalize on 
those n' flashes that end up in the past of to; then this distribution is given again by (12) 




(12) 



but with \j/ replaced by the conditional wave function 

Va^tttt; ^ M- ^^^^ 

A RELATIVISTIC GRW MODEL 

I now present a way of making the GRW model Lorentz invariant [5], exploiting the 
flash ontology and the analogy with the multi-time translation invariance. The relativistic 
model is nonlocal, excludes superluminal signalling, and converges to the GRW model 
in the non-relativistic limit c ^ o°. The version I shall present treats distinguishable 
quantum- mechanical Dirac (relativistic spin- 5) particles without interaction in either flat 
or curved relativistic background space-time. 

Let me guide you through the equations defining the model. Each of the flashes wears 
a label i E {1, . . . , A'^} corresponding to one of the A'^ variables of the wave function. The 
distribution of the flashes is again of the form (12), but based on different K operators, 
constructed in a Lorentz-invariant way. For two spacelike 3-surfaces Ei,E2, let 

[/g:L2(Zi)^L2(L2) (14) 

denote the unitary evolution (according to the Dirac equation) from Ei to E2. For every 
spacelike 3-surface L, let Az{x) be the multiplication operator on L^(S) by a Gaussian 
centered at x, 

AzWvW = f exp(-!MC^:Z))„W (15) 

where y eH, "s-dist^" denotes the spacelike distance along E, and <7 and T are the same 
constants as before. The normalizing factor ^ is chosen so that 

/ d^xA^ix) = - (16) 

Jl. T 

and may depend on y if space-time is curved. In the relativistic formula that will replace 
(5), we will use the expression A-z{xy/'^ to replace the nonrelativistic operator A(r)^/^, 
and we take E to be the hyperboloid based at the previous flash, as in Fig. 3. This 
construction cannot be done with continuous collapse as in CSL models because it is 
the (previous) flash that defines the relevant 3-surface. 

In detail, we write, for any sequence / = (jcq, JCi ,JC2, . . . , jc„) of space-time points, 

Kif) = K^„_, (xn) ■ K,, (X2) K,, (xi ) . (17) 

The collapse operator ^y(x), which we define presently, corresponds to a flash at x 
given that the previous flash was at x' , and is the analog of the nonrelativistic expression 
K^t {x) = A(r) W'f_f/, except that we use the Heisenberg picture to define it on L^(Eo), 
with Eo the 3-surface on which the initial wave function is specified (more about it at 
the beginning of the next section). We define 

KA^) = W(x')e-'-^-^'/''C/?A£W^/'f/i, (18) 



FIGURE 3. The 3-surface l,{x',x) of constant timelike distance from x' containing x 



where ^(x') means the future of V (the interior of the future Ught cone of V), |x — V | the 
timelike distance between x and x', and E = E(/,x) = {y E ^{xf) : |y — x'j — \x — x'\} 
denotes the 3-surface of constant timelike distance from x' containing x, as in Fig. 3. The 
first factor lxe-^{x') ensures that flashes for the same particle are timelike separated. The 
second factor q-\x-x'\/2z corresponds to the exponential distribution (with expectation 
t) of the proper waiting time |x — yj. The last three terms are the Heisenberg-evolved 
version of A£(x)^/^. A key property of (18) is that 



/ 



d'^xK;,{x)K^{x) = l. (19) 



Since the distribution of the next flash depends on the previous flash, we need to 
specify a seed flash xq as part of the initial condition; in fact, one seed flash Xifi for every 
particle label i. For the initial condition of the universe, they could be chosen to be, e.g., 
the Big Bang. With the notation = (xj^, . . . ,Xi^ni) and 

^f = Ilfl^Sk, (20) 

i=lk=l 

the formula for the joint distribution of the first flashes of type /, the analog of (12), 
reads 

N 



F(Xi^k e d\k ■.i<N,k< n,) = \\(S}K{fi) V 

i=l 



'd/. (21) 



By (19), this defines, at least formally, a consistent family of probability distributions. 

This completes the definition of the model. Its Lorentz invariance is manifest from 
the fact that we did not mention a preferred slicing of space-time, or any other sort of 
additional space-time structure. 



CONDITIONAL WAVE FUNCTION 



Do the seed flashes have to lie on the 3 -surface Lq on which we specify the initial wave 
function i//^? No. The measure P defined by (21) does not change if we replace Eq by any 

other spacelike (Cauchy) 3-surface Lg y/hy y/ = U^^ y/, using the unitary evolution. 
Thus, as far as (21) is concerned, Eq may lie in the past of some seed flashes and in the 
future of others, or even in the future of (some of) the random flashes governed by (21). 

This should not be understood, however, as answering the question what function is 
best regarded as the collapsed wave function on a given 3-surface £. This is the question 
I address in this rather technical section. Before, however, I will provide a couple of facts 
that help us deal with the subtleties involved in conditioning on the information that a 
certain point x was the last flash before a certain 3-surface L. First, the formula for the 
probability that no flashes occur up to the 3-surface L, i.e., the relativistic analog of (11), 
reads 



N 



(8)w.,o(r)v^ 

!=1 



(22) 



with the abbreviation 



1/2 



(23) 



where ^{T-) — l^xe'L^i^) is the future of E. And second, as a consequence of (22), 
the probability of a particular sequence fi — (x,^0:-^j,b • • • ;^i,n, ) of flashes of type / with 
Xi^i e =^(L), given that the seed flash x/^o was the last before L, is not (21) but instead 



\\®im)w\[ 

|<8)/Wx,o(S)V'| 



rd/. 



(24) 



Now I return to the question, what is the collapsing wave function (or conditional 
wave function) in this model? I have two alternative proposals. My first proposed wave 
function associated with any spacelike 3-surface L: 



(8),-w^(£)jj:(/;0y 



(25) 



where // is the finite sequence of all flashes of type i up to E, and x\ is their last one. The 
distribution of the flashes fi after Z, conditional on the flashes // up to Z, is given by 



KcmA^y'wi. 



(26) 



That is, if we want to express the distribution of the future in terms of the present wave 
function, rather than of the initial wave function of the universe, then (26) is the more 



relevant formula than (21). As we push Z to the future, i/A^ collapses whenever L crosses 
a flash, and evolves deterministically in between; the evolution law reads 



^,W,,(t)K{fi)WjZ)-WL 
¥f = 71 = ■ ^ rr , (27) 

||(g),w^,(L)i^:(/;-)w^(L)->E||' 

if E lies in the future of E, fi are the flashes between E and E, and is their last one. It 
is perhaps worth emphasizing that, for any 3-surface E, to specify Yz and the last flash 
before E (of every type) suffices for determining the conditional distribution of all events 
in the future of E given the past of E. 

My second, inequivalent definition of a collapsing wave function is 

which is not normalized, W^zW 7^ 1, and that is a feature one would not expect of a wave 
function; however, this definition leads to simpler formulas replacing (26) and (27), 
namely 



'd/ (29) 



and 



0~ = " (30) 

\\^iW,,{t)K{fi)ul'<k\\' 

The existence of inequivalent definitions of the wave function on E underlines that the 
theory is not about wave functions, but that the wave function is a tool for formulating 
the law of flashes; see [9] for further discussion of the status of the wave function. 

If we associate with every E a random wave function, be it by (25) or (28), then to this 
ensemble of wave functions there corresponds a density matrix pz. I have not succeeded 
so far in deriving a master equation for p^, and the difficulties that arise suggest that 
none exists because the evolution of p£ is not autonomous (i.e., there can be 3-surfaces 
Ei,E2 and initial wave functions i^, \]/' such that p^j = p^^ but p^^ 7^ Pzj)- ^^^^ "^^y 
part of the reason why previous work on relativistic collapse theories did not arrive at 
the present model: Diosi, for example, told me he expected that the first step towards a 
relativistic collapse model should be to find a relativistic master equation, but came to 
the conclusion that no such equation exists. 



REMARKS 

An object moving at a speed v close to the speed of light c experiences fewer flashes, 
in fact by a factor of 1 / y^l — v^/c^, as one would expect from time dilation. No 
synchronization between the flashes for different particles is postulated or obtained in 
this model. With the matter density ontology, the model would not be Lorentz invariant. 



Instead, we would need a foliation (see Fig. 2) to be able to write down a formula 
analogous to (1). 

The reason why superluminal signalling is impossible is that, by repeated application 
of (19), the marginal distribution of type-1 flashes, 

J d/2 ®K{f2) w\f = \\{K{h) ® 1) V^ll' = \x{pf''K*{f,)K{h)) , (31) 

does not depend on either H2 or 1// except through the reduced density matrix pf^ = 
tr2|v^)(VA^|. Thus, for two separated systems, the distribution of the events in one region 
of space does not depend on the external fields, which we may imagine as arranged by 
an experimenter at will, in distant regions of space, and not either on features of the 
entangled wave function beyond those captured in the reduced density matrix. 

The model is nonlocal. The basic reason is that distant flashes can be correlated, as 

||Wi)®^(/2)v^||V||Wi)V^i||'||^(/2)V^2||' unless 1^=1^1 01^2- (32) 

While not all correlations are nonlocal correlations, there is nothing in the equations of 
the model that would preclude nonlocal correlations, for example between the space- 
time locations of two spacelike separated flashes. In this model we have, rather than 
merely the nonlocal probability law for results of experiments that we have in ordinary 
quantum mechanics, a nonlocal probability law for microscopic events, the flashes. An 
interesting feature of this model's way of reconciling nonlocality with relativity is that 
the superluminal influences between correlated events do not have a direction; in other 
words, it is not defined which of two events influenced the other. This is different in 
Bohmian mechanics with a preferred slicing, where, in the case of two experiments on 
entangled particles at spacelike separation, one is earlier than the other (in the sense 
defined by the slicing) — a fact that defines a direction of the nonlocal influence. 

CONCLUSIONS 

To conclude, the model describes a possible many-particle world in which outcomes of 

experiments have (to a sufficient degree of accuracy for all cases presently testable) the 
probabilities prescribed by quantum theory. It is fully compatible with relativity, in that 
it does not rely on a preferred slicing (foliation) of space-time. While the present model 
works on Minkowski space-time and Lorentzian manifolds, a relativistic collapse model 
on a lattice (Z^ instead of M^) was developed by Dowker and Henson [23]. Another 
model, which reproduces the quantum mechanical probabilities exactly, is available 
in the framework of Bohmian mechanics [21], but uses a preferred slicing of space- 
time. Thus, with the presently available models we have the alternative: Either the 
conventional understanding of relativity is not right, or quantum mechanics is not exact. 
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